Degradation of a quantum reference frame 
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We investigate the degradation ol relerence Irames, treated as dynamical quantum systems, and 
quantify their longevity as a resource for performing tasks in quantum information processing. We 
adopt an operational measure of a reference frame's longevity, namely, the number of measurements 
that can be made against it with a certain error tolerance. We investigate two distinct types of 
reference frame: a reference direction, realized by a spin-j system, and a phase reference, realized 
by an oscillator mode with bounded energy. For both cases, we show that our measure of longevity 
increases quadratically with the size of the reference system and is therefore non-additive. For in- 
stance, the number of measurements that a directional reference frame consisting of A'^ parallel spins 
can be put to use scales as N^. Our results quantify the extent to which microscopic or mesoscopic 
reference frames may be used for repeated, high-precision measurements, without needing to be reset 
- a question that is important for some implementations of quantum computing. We illustrate our 
results using the proposed single-spin measurement scheme of magnetic resonance force microscopy. 



I. INTRODUCTION 

In quantum measurement theory, the apparatus is gen- 
erally treated as a classical system - one that is not de- 
scribed by the quantum formalism. The same is true for 
the macroscopic systems that serve as reference frames, 
for instance, a ruler with respect to which positions are 
defined, a set of gyroscopes with respect to which orienta- 
tions are defined, or a clock with respect to which phases 
are defined. Just as attempts to quantize the measure- 
ment apparatus have led many researchers to founda- 
tional puzzles (such as the quantum measurement prob- 
lem) , the quantization of reference frames has also gener- 
ated its fair share of confusion and controversy (Ref. 1] 
provides a synopsis and numerous references). 

In addition to providing, through such puzzles, an op- 
portunity for us to refine our understanding of quantum 
theory, the quantization of reference frames is also use- 
ful for answering certain practical questions. In many 
quantum experiments, one can understand certain sys- 
tems as constituting mesoscopic or even microscopic ref- 
erence frames to which other systems are compared. In 
such cases, the conventional approach wherein reference 
frames suffer no back-action may yield a poor approxi- 
mation to a full quantum treatment. 

A recent example is the investigation of how well a 
quantum optical field prepared in a coherent state can 
serve as a local oscillator for homodyne detection 0. It 
has also been proposed that the finite size of a reference 
frame can lead to effective decoherence on the system it 
describes in contexts ranging from collapse theories ^ 
to the evaporation of black holes In particular, there 
have been several investigations into the consequences of 
the quantum nature of laser fields for the manipulation_of 
physical qubits in quantum information processing \~ 



Our goal in this paper is to go beyond consideration 
of how the finite size of the reference frame affects the 
purity or coherence of systems described with respect to 
it, and to examine how this finiteness affects its longevity 
i.e., how many times it can be used in a measurement as 
an accurate reference frame. In particular, we demon- 
strate quantitatively how, if one makes measurements 
of the "orientation" of a large number of systems rela- 
tive to a single quantum reference frame (RF), then as 
a result of these measurements the state of the RF be- 
comes more mixed, more symmetric and is thereafter less 
useful for implementing the sorts of tasks for which an 
RF is required. This degradation (and ultimate deple- 
tion) is yet another reason (in addition to those provided 
in 0, H, 01 E3l) for considering RFs as an information- 
theoretic resource. 

At first glance, one might expect that the longevity of 
an RF, quantified by the number of estimations of the 
relative orientation of a system to the RF that one can 
achieve with a certain error bound before it degrades be- 
yond use, would scale linearly with the size of the RF. For 
example, if two identical directional RFs were used to- 
gether, one would naively expect the resulting combined 
RF to last twice as long as either constituent RF would 
separately. However, our analysis shows that the scaling 
does not match this expectation. In fact, the longevity of 
such a directional RF scales quadratically with the num- 
ber of spins N. We show that the same scaling holds for a 
phase reference, such as a coherent state. Specifically, the 
longevity of such a phase reference scales quadratically 
in the average excitation number (e.g., photon number). 
This result is encouraging for the potential usefulness of 
mesoscopic reference frames in quantum information pro- 
cessing applications. 

We illustrate this point with a simple calculation for 
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single-spin measurements using magnetic resonance force 
microscopy, demonstrating both that this measurement 
scheme may exhibit effects due to the degradation of its 
directional RF with repeated use, but also that it may 
be possible to use such a device for the large number of 
measurement required for quantum computing applica- 
tions without significant degradation. 

We emphasize the distinction between the quality of 
an RF, which describes how well a quantum RF approx- 
imates an ideal classical RF for the purposes of measure- 
ment, and the longevity of an RF, which describes how 
many times it can be used in a measurement whilst main- 
taining a certain quality. Whereas general uncertainty- 
principle-based arguments can often be used to deter- 
mine the scaling of the quality of a quantum RF, we 
are unaware of any such argument that can predict the 
quadratic scaling of the longevity that we determine in 
this paper. 



II. PRELIMINARIES 

We define the longevity of a reference frame as the 
number of times it can be used to perform a particu- 
lar operational task with some chosen finite degree of 
success. To explain the particular task we choose, it is 
convenient to be able to compare our RF with another, 
much larger "background" RF. We denote the RF whose 
longevity we are studying by R, which is initially cor- 
related (aligned) with the background RF. (In Ref. 
a RF R with this feature was referred to as implicated.) 
For example, R may be a gyroscope used in a laboratory 
experiment, to which directional systems are compared, 
and the background RF could be the frame defined by 
the earth or the fixed stars. 

The task to which R will be put to use is the estimation 
of the direction of a system S relative to the background 
RF. (We here use the term "direction" in a generic way 
to mean the group element relating S to the background 
RF even if the the group of transformations in which 
we are interested is not the rotation group). Such an 
estimation is achieved by measuring the relation between 
S and R, then combining the outcome with one's prior 
information about the directionality of R relative to the 
background RF, to deduce something about the direction 
of S relative to the background RF. However, because 
i? is a finite quantum system, the measurement of the 
relation between S and R causes a disturbance to the 
quantum state of R, so that after this measurement, one's 
ability to infer the directionality of another system S' is 
decreased. Thus, R can only be used a finite number 
of times for such a task before its directionality relative 
to the background RF is so poorly defined that it is no 
longer useful for this purpose. It is this degradation that 
we investigate. 

Because neither S nor R physically interacts with the 
background RF during the task of interest, we can take 
the background RF to be non-dynamical. (More pre- 



cisely, because the dynamics induced between R and S 
is completely relational, we will see that the use of the 
background RF is equivalent to simply choosing a con- 
venient gauge and has no bearing on the physics). The 
RF R, on the other hand, must be treated dynamically. 
Consequently, we assign a Hilbert space and quantum 
states to R but not to the background RF. Adopting the 
terminology of Ref. 1], R is treated as an internal RF, 
while the background RF is treated as an external RF. 
We shall also say that we are treating i? as a quantum 
RF, and the background RF as a classical RF. 

A few aspects of our direction-estimation task must be 
made specific. First, there is the nature of the measure- 
ment that determines the relation between S and R. To 
be conservative, we assume: (1) the measurement is the 
one that maximizes one's ability to estimate the relation 
between S and i?, and thus also maximizes one's ability 
to estimate the direction of S relative to the background 
RF (the figure of merit for estimation will be specified 
later) and (2) the measurement is implemented in a man- 
ner that leads to the smallest possible degradation of R 
while still satisfying feature (1). Feature (2) ensures that 
we are determining the best possible longevity for a given 
degree of success in direction-estimation. Taking into ac- 
count the fact that the more information we gain in the 
measurement, the more disturbance we create, feature 
(1) ensures that the longevity we derive will be a lower 
bound for the number of uses to which R can be put for 
any other sort of direction estimation task on S. 

The second aspect of our estimation task that must 
be made specific is the initial quantum state of R. We 
choose to investigate the state that leads to the least 
initial error in estimating the directionality of S relative 
to the background RF. We emphasize that, although the 
optimal RF states and optimal measurements may not 
necessarily be easy to achieve in practice, they establish 
the quantum limit and therefore bound any RF longevity. 

We note that an equivalent definition of longevity can 
be achieved without making any reference to a back- 
ground RF. In this case, our description of S and R would 
be subject to an effective superselection rule for the group 
that is associated with the RF For instance, if the 
RF is for orientation, then moving to a description that 
makes no reference to the background RF would imply 
adopting an effective non-Abelian superselection rule for 
the group of rotations SU(2). Although the usefulness of 
this mode of description has been emphasized elsewhere 
[Hi, we shall not make use of it in the present paper. 
In what follows, we assign quantum states to S and R 
that are non-invariant under the group of interest, even 
though the measurement of the relation of S* to i? is in- 
variant under this group. That is, the quantum states are 
assigned relative to an arbitrary background RF, which is 
not used in the measurement. One may liken this proce- 
dure to choosing a gauge for convenience when describing 
gauge-invariant processes. 
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III. DIRECTIONAL REFERENCE FRAME 

In this section, we study the degradation of a direc- 
tional RF, that is, a reference for a direction in space. 
(Note that a directional RF does not provide a full Carte- 
sian frame, as it does not provide a reference for rotations 
about its axis.) 

A. Measurement of a spin-1/2 system relative to a 
spin-j directional RF 

We use a spin-j system for our quantum directional 
RF, with Hilbert space Mj. The initial quantum state of 
the spin-j system is denoted p'-*'-' . (Later, in Sec. IIII CI 
we will determine the state that minimizes the error in 
the estimation of direction.) We describe this RF relative 
to a background RF, as noted in Sec^ and choose it to 
be aligned in the +z direction relative to the background 
RF; we emphasise, though, that this alignment with a 
background RF is essentially a choice of gauge. Because 
the quantum RF will serve only as a reference direction, 
and not a full frame, we can choose it to be invariant 
under rotations about the z axis without loss of general- 
ity. Thus, the initial quantum state, p*^*^^ is diagonal in 
the \j,m) basis of Hj, that is, the basis of simultaneous 
eigenstates of and Jz ■ 

The systems to be measured against the quantum RF 
will be spin-1/2 systems, each with a Hilbert space Hi/2- 
We choose the initial state of each such system to be the 
completely mixed state 1/2, where / is the identity op- 
erator on M.i/2- This corresponds to maximal ignorance 
about the spin-1/2 system. Our quantum RF will be used 
to measure many such independent spin-1/2 systems se- 
quentially. We shall assume trivial dynamics between 
measurements, and thus our time index will simply be an 
integer specifying the number of measurements that have 
taken place. The state of the RF following the nth mea- 
surement is denoted p^"-*, with p*^*^^ denoting the initial 
state of the RF prior to any measurement. We consider 
the state of the RF from the perspective of someone who 
has not kept a record of the outcome of previous measure- 
ments. Thus, at every measurement, we average over the 
possible outcomes with their respective weights to obtain 
the final density operator. 

A measurement of the relative orientation of a spin-1/2 
particle to a spin-j system is represented by operators 
that are invariant under collective rotations. The mea- 
surement that provides the maximum information gain 
about the relative orientation between a spin-j and a 
spin-j' system has been determined in 12]. It is sim- 
ply a measurement of the magnitude of the total angular 
momentum J^. That is, it is the projective measure- 
ment {Hj}, where Hj projects Mj^Mj' onto the total 
eigenspace with eigenvalue J{J +1). For the case where 
j' = 1/2, the optimal measurement for determining the 
relative orientation is represented by the two-outcome 
projective measurement {11+ = nj+i/2,H_ = Tlj-1/2} 




B. Measurement-induced update of the directional 
RF state 

A given projective measurement can be associated with 
many different update maps. As discussed in Sec.^ we 
choose the update map that is minimally disturbing. In 
App.0 it is shown that this corresponds to adopting the 
Liiders rule for updating: for a measurement outcome 
associated with the projector H, the quantum state p is 
updated to npH/Tr(pH). 

Thus, under the optimal measurement for relative ori- 
entation, the evolution of the quantum RF as a result of 
the nth measurement is 

=£(p(»)) (1) 

where 

£(p)-Trs( n,(p®//2)H,) , (2) 

ce{+,-} 

and Trs denotes the partial trace over the spin-1/2 sys- 
tem. 

The map £ can be written using the operator-sum rep- 
resentation 13] as 

^^p) = \ E E ^>^aL (3) 

ce{+,-} a,fce{oa} 

where 

E:,^{a\n,\b) (4) 

is a Kraus operator on Hj and {|0) , |1)} is a basis for 
These operators can be straightforwardly deter- 
mined in terms of Clebsch-Gordon coefficients. 

We note that the evolution of the RF is not unitary 
(and would not be unitary even if the measurement re- 
sult was kept), because the quantum RF becomes en- 
tangled with the measured system which is subsequently 
discarded. Thus, the map £ is irreversible and conse- 
quently one may say that the quantum state of the RF 
undergoes decoherence. 

In App. ^ we provide a recurrence relation for the 
matrix elements of the state p'-"-' . 

C. Measure of directional RF quality 

One simple measure of how much the RF degrades is 
the fidelity 13j between the state after n measurements, 
p'^"\ and the initial state, p'^'^\ which clearly decreases 
with n. Another natural measure is the asymmetry of 
the state which also is found to decrease with n. 
However, rather than use these measures, we will instead 
choose to quantify the quality of the RF operationally as 
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the average probability of a successful estimation of the 
orientation of a spin-1/2 system in a pure state (relative 
to the background RF). For simplicity, we assume that 
this fictional "test" spin-1/2 system is with equal proba- 
bility cither aligned or anti-aligned with the background 
RF. It is worth making special note of the fact that in 
the sequence of measurements that cause the RF degra- 
dation, the state of the spin-1/2 system is not assumed 
to be a pure state that is aligned or anti-aligned with 
the RF. It is assumed to be completely mixed. However, 
after each measurement we ask: if a spin-1/2 system in 
such a pure state were compared to the RF, how well 
could we estimate whether it was aligned or antialigned? 

Denote the pure state of the test spin-1/2 system that 
is aligned (anti-aligned) with the initial RF by |0) (|1)). 
For a spin-1/2 system prepared in the state |0), the prob- 
ability of success, that is, of finding the correct (j + 1/2) 
result, is 

F,"(n) = Tr([p(") ® |0)(0|]n+) = TTnip^^^^ E+,) . (5) 

Similarly, for a spin-1/2 system prepared in the state |1), 
the probability of the correct (j — 1/2) result is P^in) = 
Trij(p(")£;fi). Thus, the average probability of success is 

P,{n) = ^TTn{p^-HE++E^,)). (6) 

Using Clebsch-Gordon coefficients to determine the 
Kraus operators of Eq. , we find 

j 

^s(n) = l-(5^ J2 (j-m+l/2)Trj^(p(")|j,m)(j,m|). 

m=-j 

(7) 

We remind the reader that, at every measurement, the 
resulting density matrix for the RF has been averaged 
over the possible previous outcomes with their respective 
weights. Thus, this measure Ps{n) is applicable to situ- 
ations where the prior measurement record is discarded. 
In addition, this method also quantifies the average case 
longevity even if the measurement record is kept. 

We now focus our attention on a particular initial state 
of the RF. We choose the state that is optimal for our 
measure of RF quality, in other words, we choose the p'-*'-' 
that maximizes -Ps(O). This involves minimizing the sum 
in Eq. (jT)), which may be written as 

Tr(iVV)), (8) 

where 

j 

M= J2 {j-m + 1/2) \j, m) (j, m| . (9) 

Clearly then, we should choose to be an eigenvector 
of M associated with the smallest eigenvalue, namely, 
= m = j) {j, m = j \ . Thus we see that the quan- 
tum state that maximizes our measure of quality of 



the RF is an SU(2) coherent state, yielding Ps(0) = 
l-l/[2(2j + l)]. 

(This result is not in conflict with that of Peres and 
Scudo |15| . where it was demonstrated that the prod- 
uct state of N parallel spins (which is equivalent to the 
SU(2) coherent state) was not the optimal state for the 
transmission of a reference direction, because their sys- 
tem was not assumed to be a spin-j representation, as we 
have done here. An interesting direction for future work 
is to extend the results of in order to determine the 
optimal measurements of spin-1/2 systems against these 
more general RF systems and to explore the degradation 
in this case.) 

D. Scaling of the longevity of a directional RF 
with size 

Our problem is to solve Eq. |^ to determine the RF 
state as a function of rt. In the limit j 3> 1, one can 
obtain a closed-form expression for the diagonal elements 
of the state of the RF at step n, assuming an initial state 
p'"^ = It is 

Tr(p(")|j, m)(j, m|) 

= (-1)^-" f . " ) U-my.i-2j)U{j-m+l, n+2, -2j) , 

(10) 

where U{a, b, z) is Rummer's U function Details are 
given in App. IbI 

Using this result, we find the rate at which Pa{n) de- 
creases with n to be 

?^(")-^+2^(l-W^)"- (11) 

We plot this for various values of j in Fig. ^ 

The initial slope of this function bounds the rate of 
degradation. It is 

R = Ps{l) - Ps{0) = -2j/{2j + if . (12) 

Thus, in the large j limit, we have the rate of degradation 
with n satisfying R > — l/(4j^). 

Let e < 1 be a fixed allowed error probability for the 
spin-1/2 direction estimation problem. After n measure- 
ments, the probability of successful estimation is lower 
bounded by 1 -|- nR, so the number of measurements re- 
quired to ensure that this bound be greater than 1 — e is 
— e/i?. Consequently, the number of measurements that 
can be implemented relative to the spin-j RF with prob- 
ability of error less than e is 

"max ^ . (13) 

This result implies that the number of measurements for 
which an RF is useful, that is, the longevity of a RF, 
increases quadratically rather than linearly with the size 
of the RF. 
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FIG. 1: A plot of the success probability for a test particle to 
be measured correctly against the directional reference frame, 
versus the number of previous spins that have been measured 
against the frame. The plots (from bottom left to top right) 
are for a reference frame of spin j = 1, 3, 5, 10, 15, 25, 50, 100, 
respectively. 



FIG. 2: A loglog plot of the longevity (number of measure- 
ments before the error probability for a test particle exceeds 
a fixed threshold) of a directional reference frame versus the 
size of the frame. Plots for various different error thresholds 
are given. A clear quadratic asymptotic dependence emerges 
despite the relatively small reference frames indicated: linear 
regression yields gradients which all lie within 2 ± 0.1. 



Numerical calculation of the longevity's dependence on 
RF size for various choices of error threshold confirm this 
result; see Fig. El 

This result imphes that, in order to maximize the num- 
ber of measurements that can be achieved with a given 
error thereshold, one should combine all of one's direc- 
tional RF resources into a single large RF and perform 
all measurements relative to it. It is also worth noting 
that if the directional RF is composed out of N spin- 1/2 
systems, so that it has size j — N/2 + 1, the number 
of measurements that can be performed with probabil- 
ity of error less than e scales as iV^ , whereas one might 
have naively thought that it would be good for only as 
many measurements as there are constituent spins. Fi- 
nally, note that an SU(2) coherent state for N spin-1/2 
systems is a product state so that entanglement does not 
appear to be responsible for this property. 



E. Application to single-spin measurement using 
magnetic resonance force microscopy for solid-state 
quantum computing 

The physical nature of measurement apparatuses and 
reference systems is particularly relevant for quantum 
computation. The stringent requirements placed on mea- 
surement devices for use in quantum computing include 
(i) large coupling strengths, often requiring the measure- 
ment device to be made very small and placed close to 
the quantum registers (qubits); and (ii) the apparatuses 
must be well-described by classical, noiseless devices in 



order to obtain the required levels of accuracy (e.g., per- 
form projective von Neumann measurements). These two 
requirements appear at first glance to be mutually exclu- 
sive, but may be satisfied by novel mesoscopic measure- 
ment schemes. Consider, as an example, the proposal 
to use magnetic resonance force micro scopy ( MRFM) to 
measure the direction of a single spin [isl Il9l | in a solid- 
state quantum computer. This scheme makes use of a 
small magnet placed on the tip of a nano-mechanical res- 
onator; monitoring the oscillations of the resonator can 
be used to measure the direction of a spin placed near 
the magnet. To couple strongly to a single spin, the 
measurement device must be cold and very small (i.e., a 
resonator with a mass of the order of picograms, with a 
magnet on the tip of length scale 10 nm [l^)- However, 
objects on this scale are poorly approximated as classi- 
cal, macroscopic measurements apparatuses, and thus a 
semi-classical or full quantum treatment of the measuring 
device is required. 

The results we have presented here can be applied to 
the proposed use of MRFM to measure the direction of a 
single spin, and its associated applications to solid-state 
quantum computing. The design parameters of a scheme 
to obtain single-spin measurement are outlined in |l9j| . 
in which the magnet consists of an iron sphere of diam- 
eter ~ 300 A. Treated as a quantum system, then, the 
magnet can be approximately described by 10^ paral- 
lel spins. With an error tolerance of 10"'* (a reasonable 
threshold for quantum computation) this RF would per- 
form 10^ optimal measurements before degrading beyond 



6 



further usefulness. Thus, if the measurement scheme can 
be chosen such that the disturbance on the RF due to 
each measurement is comparable with that of the opti- 
mal update map, it appears that a mesoscale RF such 
as this may be suitable for repeated use in large-scale 
quantum computation. 

IV. PHASE REFERENCE 

We now consider the case of a phase reference. Our 
systems will consist of oscillator modes, taken to be opti- 
cal modes for clarity. It should be noted that we focus on 
the optical case only to make the description of the re- 
sults less abstract; the results are applicable to any phase 
reference. 



A. Measurement of single-rail qubit relative to 
phase reference 

For our phase reference, we take a single oscillator 
mode with Hilbert space denoted by M^. (Note that 
we could instead use a multimode phase reference, or 
even many qubits prepared as "refbits" Again, the 

initial quantum state of the RF is denoted p^^'' . We de- 
scribe this RF relative to a background phase reference, 
and arbitrarily choose its phase to be zero (our choice of 
gauge). 

The systems to be measured against the RF will be 
single-rail qubits, i.e., oscillators with states restricted 
to the lowest-energy 2-dimensional subspace spanned by 
Fock states |0) and The Hilbert space is denoted H5. 

We consider the particular estimation task wherein the 
system is promised to encode phase or phase tt with 
equal probability, which is to say that it is promised to be 
in a state \+) or |-), where |±) = 2'^/^ (|0) ± |1)), with 
equal probability. In this case, an optimal measurement 
for estimating whether the relative phase between system 
and quantum RF is or tt is the two-outcome projective 
measurement {11+, n_} with 

00 

n+ = ^ |m, +> (m, +\ + \0)r (0| ® \0)s (0| 

m— 1 

00 

n_ = ^|m,-)(m,-|, (14) 

m— 1 

where 

|m,±) = -l={\m)n\0}s ± |m-l)H|l)s) ■ (15) 

This is demonstrated in App. [21 

It is worth noting that this measurement is a coarse- 
graining of the optimal projective measurement of 
phase on each of the 2-dimensional eigenspaces of 
total photon number; we further note that this measure- 
ment was used in 21] for the same purpose. 



We will be interested in determining how the RF state 
evolves after many such measurements on distinct sys- 
tems. For simplicity, we consider the RF state as it is 
described by an observer who does not know whether the 
system was in the state |-|-) or |— ) initially, and because 
the two are presumed equally likely, the initial state of 
the system is presumed to be the completely mixed state 
1/2. We also assume that the observer does not keep 
track of the results of the measurements. 



B. Measurement-induced update of the phase 
reference 

As with the directional example, the update rule 
for this two-outcome projective measurement that is 
minimally-disturbing is the Liiders rule, described in 
App.m 

With this choice of update rule, the evolution of the 
quantum RF as a result of the nth measurement has the 
same form as Eqs. (|1I4|) but where the E^j^ operators are 
now given by 

E+^h,{lR + \0){Q\), E-, = ^{In-\0){0\), (16) 
E+=E-, = ^In, (17) 

E+ = -i?r„ = {E+y = -iE,,y = , (is) 

where Ir is the identity on M^^ and 

00 

A=Y^ \m){m + l\. (19) 
The update map due to a single measurement is therefore 

S{P) = ip+ i|0)(0|p|0)(0| + \A^pA+lApA^ . (20) 

Again we note that although it may be difficult to 
achieve such an update map in practice, it is nonetheless 
interesting to consider because it defines the quantum 
limit of RF longevity. If the systems under consideration 
are modes of the electromagnetic field, this update map 
is particularly impractical given current technology be- 
cause it corresponds to a quantum non-demolition mea- 
surement on the electromagnetic field. An interesting 
problem for future research is to determine how much 
degradation, in excess of the quantum limit, is incurred 
by the use of more realistic update maps. 

C. Measure of phase reference quality 

As with the directional case, we operationally quantify 
the quality of the RF at step n as the average probabil- 
ity of success in a hypothetical measurement of whether 
a qubit system is in phase or out of phase with the back- 
ground phase reference, that is, whether the qubit is in 
the state or |— ) , given equal prior probability for 
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these two possibiUties. Denoting the state of the RF at 
step n by p'-"-', and using the optimal measurement (|14|l . 
the average success rate is found to be 



1 



2 ^ 

m=0 



Re 



/ (n 



m+1 



(21) 



where = Tr(p(")|m)(m + 1|). 

Just as we did for the directional case, we consider the 
degradation in the case of an initial state of the RF that 
is initially optimal for our measure of RF quality, i.e., the 
state that maximizes Ps(0). However, if the accessi- 
ble Hilbert space is taken to be infinite-dimensional, 
then this probability can trivially be made arbitrarily 
close to one. To obtain a non-trivial and physical result, 
we must bound the energy of the state in some way; 
we may then use this bound as a measure of the "size" 
of the RF, analogous to the use of total spin j for the 
directional case. 

One natural choice for implementing this bound is to 
limit the maximum photon number to some fixed value 
N, i.e., demand that p'-"-' has support entirely on the 
(N + l)-dimensional subspace spanned by {\m) , m = 
0, . . . , N}. We now seek the state p'^''^ that maximizes 
Ps(0) subject to this constraint. This requires us to max- 
imize the expression 



Tr(AW"^) 



(22) 



where 



N 



Ma,= ^(|m)(m + l| + |m + l)(m|). (23) 



Thus we must choose p*-"-* to be the eigenvector of Mm 
with minimal eigenvalue. The solution is 



N, 



N 

= A/" sin 



(m -I- l)7r 



|m) , 



(24) 



where A/" is a normalization factor. This is the same 
state that is optimal for determining the relative path 
length between two arms in an interferometer, for several 
different figures of merit ^23\ . We provide a derivation of 
Eq. H24|l in App. ^ For this state, we have 



Ps{0) = -{1 + cos{tt/{N + 2))) 



1 



2N^ 



for large N . 



(25) 



We note, however, that our method of bounding the 
energy of the initial state p^°^ excludes states, such as 

the coherent state \a) = e^l"l X]^^o('^"/"^^)l"')' 
have non-zero support on the entire Hilbert space Mfj but 
nevertheless have a finite energy. It is beyond the scope 
of this paper to explore all the varied possibilities for de- 
termining optimal phase references for alternate choices 



of bounds. However, because of the ubiquity of coher- 
ent states as phase references in many physical systems, 
we will also consider how the longevity of a coherent 
state compares with the optimal phase-encoding state 
with bounded number 1241) . We compare the size of such 
states through average photon number n, which is given 
by jap for a coherent state and by N/2 for an optimal 
phase-encoding state containing at most N photons. 



D. Degradation of phase reference 

To determine the rate at which Psin) decreases with 
n, given an initial state p*^"' for the RF, we need only 
determine the evolution of the matrix elements Pm,m+i- 
To this end, we use Eq. (^0)) to obtain the following re- 
currence relation for these elements, 



1 



1 



1 



in) 



m,m-\- 



(26) 



As it turns out, the elements on the m, m -f- 1 diagonal 
depend only on elements on this diagonal at the previous 
time step. The solution to this recurrence relation is 
found to be 



m,m+l 



n+m+1 

■ E ' 



\( 2n \ _ ( 2n \ 



This expression can be used in Eq. (|21|l to determine nu- 
merically how Ps(n) decreases with n for a given choice 
of initial state p'^^'. In Fig. |31 we plot Ps{n) for an ini- 
tial optimal phase-encoding state with number bound N, 
Eq. H24|l. for various values of N. For comparison, we also 
plot Ps(n) for an initial coherent state with mean num- 
ber jap = N/2 (so that we are comparing states with the 
same mean number of photons). We note that, although 
the optimal phase-encoding state with bounded number 
gives a superior initial probability of success Ps{0) when 
compared to that of the coherent state, the latter appears 
more robust and does not degrade as quickly. This result 
demonstrates, perhaps surprisingly, that optimizing the 
initial quality of a quantum reference frame (quantified 
by the initial probability of success) does not optimize 
the longevity of the state. It is an interesting line of fu- 
ture research to determine which quantum states give the 
optimal longevity. 

We find numerically that the longevity of a phase ref- 
erence, both for the optimal phase-encoding state with 
bounded number and for the coherent state, scales as n^, 
the mean photon number squared; see Fig.^ Specifically, 
we find that Psin) reaches a fixed value independent of 
n after a number of uses equal to n^. Again, we find that 
the number of measurements for which a phase RF is use- 
ful, that is, the longevity of a RF, increases quadratically 
rather than linearly with the size of the RF. 
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FIG. 3: A plot of the success probability for a test mode to be 
measured correctly against the phase reference frame, versus 
the number of previous oscillators that have been measured 
against the frame. Plots are shown for both coherent states 
and the optimal phase-encoding states for bounded number. 
The plots (from bottom left to top right) are for a reference 
frame of average photon number m = 1,4, 9, 16, 25, 36, 49, 64, 
respectively. The fact that the plots for the coherent state 
and the optimal phase-encoding state are nearly overlapping 
at m = 4 and 9 is a result of the similarity of the number 
distributions of the two states at these values of m. 




1 1.5 2 2.5 3 3.5 4 



log( mean photon number ) 



FIG. 4: A loglog plot of the longevity (number of measure- 
ments before the error probability for a test particle exceeds 
a fixed threshold) of a phase reference frame versus the size 
of the frame. Plots for various different error thresholds are 
given. A clear quadratic asymptotic dependence emerges de- 
spite the relatively small reference frames indicated: linear 
regression yields slope fits which all lie within 2 ± 0.1. 



V. DISCUSSION 

We have presented results on the degradation of both 
a directional and a phase reference. In each case we have 
essentially uncovered a quadratic scaling of the number 
of uses the frame can endure in terms of the size of the 
frame. An open question is whether this quadratic scal- 
ing of longevity holds for all types of RFs. 

It is interesting to consider whether a semi-classical 
model would suffice for describing the degradation of an 
RF. In particular, one could employ a description of a 
dynamical classical RF that uses a probability distribu- 
tion on the space of "orientations" which encodes ones 
ignorance of the exact direction, and also use a mea- 
surement theory that relates the probability of successful 
measurement to the uncertainty (variance) of this distri- 
bution. This uncertainty would increase with repeated 
measurements; a simple model might be to use diffu- 
sion (i.e., a classical random walk) to characterise this 
increase. Our key result in this paper is to determine the 
rate of degradation and its dependence on the size of the 
RF. If this rate was to be determined solely from a semi- 
classical model, then the model would require a sophisti- 
cated measurement theory for determining the direction 
of a quantum spin relative to a classical object of a given 
size. We believe that a heuristic, semi-classical model 
that predicts the quadratic scaling of the number of uses 
of an RF in terms of its size would be very desirable. 

It is worth noting that we have considered only the con- 
sequences of measurements on an RF, and not the conse- 
quences of a continual interaction with an environment. 
In practice, RFs do interact with their environments, and 
so the impact on their longevity is an important question 
for future research. It should be noted that although 
many environments are likely to act in a manner similar 
to measurements, thereby increasing the rate of degra- 
dation, others might in fact act to continuously realign 
one's quantum RF with another, larger, RF which exists 
in the background. For instance, if the RF is composed 
of a number of spins, and is continuously interacting with 
a strong uniform magnetic field, then thermalization will 
tend to cause the spins to become aligned with the field. 

Spatial/directional RFs are ubiquitous, and so one 
might question the need for considering their degrada- 
tion. It could be argued, for example, that any direc- 
tional RF which has degraded can be realigned with some 
background RF in the lab (treated as a macroscopic di- 
rectional RF) at any time. However, such realignment 
may prove to be difficult for microscopic apparatuses 
that are placed close to the quantum systems to be mea- 
sured. In addition, realignment of a degraded shared 
RF costs resources - an issue that can be important es- 
pecially within the context of quantum communication. 
Our main result shows, however, that such realignment 
may not be necessary for many situations, such as the 
MRFM example discussed above. 

Although some RFs, such as directional RFs, are in- 
deed ubiquitous, RFs for many other degrees of freedom 
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are not. Some RFs must be painstakingly prepared via 
some controlled quantum process. One example is the 
BCS ground state of a superconductor, which serves as a 
phase reference for experiments involving superconduct- 
ing qubits; another example is a Bose-Einstein conden- 
sate. An analysis similar to the one presented in this 
paper for any such frame could be used to determine the 
RF size necessary to keep errors below a specified thresh- 
old for a specified number of uses. 

Although in this paper we have considered the amount 
of degradation that results from measurements that max- 
imize the information gain about the relative orientation 
of system and RF, it is clear that there is a trade-off 
between longevity and information gain; the more infor- 
mative the measurement, the more degradation to the 
RF. An interesting problem for future research is to de- 
termine the precise nature of this tradeoff. 
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APPENDIX A: PROOF THAT LUDERS RULE 
UPDATING IS MINIMALLY-DISTURBING 

A measure of the quality of a RF, henceforth a measure 
of "frameness" , is operational if it quantifies the maxi- 
mum success with which one can achieve some task in a 
variation over all protocols that use only the resource and 
invariant operations. Any operational measure of frame- 
ness F is by definition a monotone with respect to oper- 
ations that are invariant under the action of the group 
associated with that RF, that is, if p is the quantum state 
of a RF for the group G, then 



F{£{p)) < F{p) . 



(Al) 



for all £ that are invariant under G, that is, for which 
£ [U{g){-)W{g)\ = U{g)£{-)W{g) for all g e G. We call 
such measures frameness monotones. 

The maximum probability of a successful estimation of 
a direction of a system relative to a quantum reference 
frame is clearly an operational measure of frameness (it 
involves a variation over all estimation protocols that use 
invariant operations and the RF resource). It is therefore 
by definition a frameness monotone. 

We now make use of the following proposition con- 
cerning CP maps that are update maps for a projective 
measurement. 



Proposition: For any map £ such that £'^{I) = 11, 
where 11 is a projector, we have 

£{p) = £{UpU) . (A2) 

Proof. Any Krauss decomposition of £ is of the form 

£{p)=Y,K^pKl, (A3) 



with 



Y,KlK,=Ii. 



(A4) 



Given that every term in the sum is positive, that is, 
KlKf, > 0, it follows that 



supp{KlK^) C supp(n) , 



(A5) 



where supp(A) denotes the support on the Hilbert space 
of the operator A. It follows that 



snMyKlKi,) C supp(n) , (A6) 

and thus 

^kIk^ = ^kIk,U. (A7) 
But now using the polar decomposition = 



U^^J KjiKi^, we infer that 



(A8) 



Our claim then follows trivially. QED. 

Now, by virtue of the fact that £ is invariant under G, 
and by virtue of the fact that the probability of successful 
estimation of a direction, Pg, is a frameness monotone, it 
follows that 



P,{£{p)) = Ps{£{IlpU)) 

< -Ps(npn) 

Normalizing our density operators, we have 

Ps(J#.)<Ps(^). 



'Tr{Ilp)' 



^TV(np)' 



(A9) 
(AlO) 



(All) 



Thus the probability of successful estimation given a 
Liiders rule collapse map p Ilpll/Ti{llp) is an up- 
per bound for the probability of successful estimation for 
any collapse map associated with 11. Thus one can do 
no better than to use the Liiders rule collapse map. 



APPENDIX B: DERIVATION OF THE 
RECURRENCE RELATION FOR A 
DIRECTIONAL RF 

In order to explicitly solve for the degraded state p'"^ 
of a directional RF as a function of n, we now derive a 
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recurrence relation for the matrix elements of the state. 
The E^f^ operators can be obtained in the |j, to) basis 
using Clebsch-Gordon coefficients by making use of the 
decomposition 



n+ = 



E 

M=-(i±i) 



|j±i,M)(j±i,M| 



(Bl) 



where 



|j±iM) 



RS 



(j, i;i,ib-±i,M)|j, M 



^>«|0)5 



(B2) 



Denoting (j, to|^ i? |j, to')^ by [i^Jxa:', where x = j — m, 
we have 



2j+l "'^'V 



(£+1), 

2j+l 
2j+l "xy 



where < x, ?/ < 2j. 

Using these expressions, we are led to the following 
recurrence relation: 

(2j + l)2p(«+i) = [2f + 2j + 1 + 2(j - - 

+ ^xy{2j + l-x){2j + l- 

+ v/(a; + + l)(2j - x){2j - y)pi%y+, . (B3) 

Note that if the initial state is diagonal in the \ j, to) basis 
(as is noted in Sec. IIII Al we can demand this without 
loss of generality), it remains diagonal. In this case, the 
recurrence relation is simply 



(2j + 1) ^("/'^ = (2j' + 2J + 1 + 2(j - X 



+ x{2j + l-x)p'-^'\^_^ 
+ {x + l)i2j-x)pi%,,^,. (B4) 

We now explicitly solve for p^") for the optimal initial 
state = for j ^ 1. The recurrence relation 

in this limit is 



^ir'^-[l-IT(2a: + l)]p; 



(") 



+ i] [^p.-i..-i + + i)pr+i,x+i] • (B5) 

Expressing the initial condition in matrix form as p^qq — 
1 and all other elements equal to 0, we find the diagonal 
matrix elements after n measurements to be 



pin) 
rxx 



k+x 



k\ k\ 



k)\x) {2jY ' 



(B6) 



which can be expressed in closed form as Eq. (|10|l . One 
can verify that this is indeed a solution to the recurrence 
relation. 



APPENDIX C: OPTIMAL MEASUREMENT OF 
RELATIVE PHASE 

We seek the optimal measurement for estimating 
whether two oscillators are in phase, or tt out of phase, 
assuming the second has a maximum occupation number 
of 1 (i.e. the second oscillator is a qubit), and assum- 
ing that the and tt relative phases occur with equal 
probability. 

The total Hilbert space that we need to take account 
of is 



0H„® (Ho 



e 



n=0 



m=0 



(CI) 



where 



{span {|to) |0) , \m- 
span{|0)|0)} 



-1>|1>} 



for TO > 
for TO = 



By Schur's lemma, any positive operator that is invariant 
under collective phase rotations must be block diagonal 
with respect to ®„i^o^»"- other words. 



Ex = 



OO 

EE 

rn— ^ 



w^^^ G 



(C2) 



where {Grn,p,}p, is a POVM on K™, and where > 



and J2> 



that Wm'!f_i is a probability distribution over A and conse- 
quently that the Ex can be obtained by random sampling 
of the Gm fi- Thus, we may as well simply measure the 
POVM {G„^^}„_^. 

Because the measurement outcome associated with 
the POVM element Kq yields no information about the 
relative phase, we adopt the arbitrary convention that 
upon obtaining such an outcome a relative phase of is 
guessed. 

Now consider the measurement outcomes associated 
with POVM elements confined to for to > 0. Within 
K„j, the optimal POVM for estimating whether the rela- 
tive phase is or TT may be assumed to be covariant with 
respect to the group of relative transformations. Because 
there are only two transformations, namely, V{Q) — I® I 
and T^(7r) = I(E)e™^ , there need only be two POVM ele- 
ments, which we denote Gm.o and Gm,Tn and covariance 
implies that 



1 to ensure that '}2,x-^>^ ~ ^- follows 



Gn 



V{TT)Grn,0V{TT)^ 



The optimal POVM may also be assumed to be rank 
1 in Km, so that 



where 



Gm,0 = \gm) (ffml • (C3) 



l5m) = gm,o \m) |0) -f g^.i |m - 1) |1) . (C4) 



11 



Noting that 

V{tt) |g,n) = 5m,o |m) |0) - g,„,i |m - 1) |1) , (C5) 
it follows that 

Gm.Q 

+ 2|g„al'k-l)(m-l|® |1)(1|. (C6) 

But given that Gm.o + G^.tt = -^k„ , we conclude that 
l5m,o|^ = \gm,i\^ = 1/2. Thus 

2-1/2 (|„i)|0)+e^*™|"^-l)|l» ^ 

for some phase (^m- Because we have assumed that the 
outcome associated with the POVM element Go is the 
one which leads to a guess of relative phase, the op- 
timal POVM must have (pm — 0. Thus, \gm) = |"^, +) 
and V{'K)\gm) = I'ti, — ), where |r7T,,±) is defined as in 
Eq. 

Coarse-graining all the POVM elements that lead to 
a guess of relative phase (which includes the projector 
onto |0) |0) according to the above convention) into a 
single POVM element 114., and all those that lead to a 
guess of TT relative phase into a POVM element n_, we 
find that the optimal POVM has the form of Eq. l(Tt| . 



where In is the identity operator on the space of N or 
fewer photons. Defining Gn = Afjv — A/at, one finds that 



detGAT = -AdetG 



detG 



JV-2 



(D2) 



for which the solution is 



detGw = L/Ar(-A/2) , 



(D3) 



where the Un are the Chebyshev polynomials of the sec- 
ond kind, given by L/,„(cos0) = • Given that 
Un{x) — ztU]\i{~x), it follows that the characteristic 
equation is C/Ar+i(A/2) = 0, and thus the largest eigen- 
value is A™^'' = 2cos(7r/(iV + 2)). 

To find the eigenvector associated with the largest 
eigenvalue, we must solve Mn\iP) = X^^^li/;). Defining 
V^m = (m|V'), we have ^m+i + ipm-i = X^'^^tpm for 1 < 
m < N — \. At m = 0, we have V'l = A^'^'^V-'o, and at 
m — N, we have ipN-i = X"^^^ipN- The solution is 



APPENDIX D: OPTIMAL RF STATE FOR 
ESTIMATING RELATIVE PHASE 



V'™ = t/m(A'"'^72)=AAsir 



(m -f l)7r 



iV + 2 



(D4) 



The characteristic equation we must solve is 

det(Mjv A/w) = 0, (Dl) 



where A/" is a constant. The coefficients i/^m fall to zero 
at TV -I- 1, verifying the presence of a cut-off in photon 
number at to = A^. This solution confirms Eq. (|24|) . 
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